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Preface 


This  report  is  the  result  of  my  attempt  to  implement  a 
method  of  detecting  the  single  visual  evoked  response  (VER). 
The  problem  is  analogous  in  many  respects  to  the  radar  de- 
tection problem.  That  is,  a small  signal  containing  the  use- 
ful information  is  buried  in  a higher  amplitude  noise  back- 
ground. There  are,  however,  a considerable  number  of  vari- 
ables in  evoked  response  detection  not  all  of  which  are  fully 
understood.  The  relationship  of  the  signal  (VER)  and  the 
noise  (on-going  electrical  activity  of  the  brain)  is  a prime 
example. 

Most  of  the  visual  evoked  response  research  to  date  has 
dealt  with  the  averaged  VER.  This  requires  the  presentation 
of  a large  number  of  stimuli  and  the  computation  of  an  aver- 
aged response.  There  are  possible  applications  of  the  VER 
that  require  positive  detection  after  a single  stimulus  pre- 
sentation. It  was  toward  this  end  that  this  project  was 
undertaken. 

I wiBh  to  acknowledge  my  indebtedness  to  Captain  Gregg  L. 
Vaughn  for  his  invaluable  help  during  this  study.  I also 
wish  to  express  my  sincere  appreciation  to  Major  Robert 
O'Donnell  for  his  assistance  during  and  after  the  data  acqui- 
sition phase.  I also  give  my  thanks  to  Jack  Capehart  for 
both  his  work  in  digitizing  the  analog  data  and  his  program- 
ming assistance.  Finally,  a special  thanks  is  accorded  to 
my  wife  and  three  sons  for  their  patience  and  understanding 
during  this  period. 
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Abstract 


^ The  visually  evoked  response  (VER)  is  one  neural  response 

(^JP*  JlOU^kj  lOltl* 

to  an  external  eventr]  Through  computer  averaging  techniques 


the  evoked  response  to  a particular  stimulus  is  readily  avail- 
able, However ,>there  are  some^anticipated  applications  which 
require  detection  of  the  VER  after  a single  stimulus  presenta- 
tion. a -D?r ook  ^ i*  "Hi*  t-p^>orT  _ ^ 

This  report  describes  an  attempt  to  detect  the  single 

N 

visual  evoked  response.  The  approach  taken^is  based  on  the 
similarity  of  the  problem  to  that  of  radar  detection  in 

which  a known  signal  is  buried  in  noise.  Using  'the  discrete 

• 

Fourier  transform  rBPTv  and  other  digital  signal  processing 
techniques.,  a filter  matched  to  the  VER  template  was  computed 
and  implemented  on  a digital  computer. 

Computation  of  the  matched  filter  required  an  estimate 
of  the  noise  power  spectrum,  where  the  noise  in  this  case  is 
the  on-going  electrical  activity  of  the  brain.  This  was  ac- 
complished by  averaging  modified  periodograms . Data  from 
both  the  vertex  and  occipital  regions  were  filtered  to  deter- 
mine if  one  region  produced  better  results  than  the  other. 


The  results  indicate  that  the  matched  filter  will  detect 
the  single  evoked  response.  However,  the  false  alarm  rate 
ia  high.^  The  filtered  output  for  data  from  both  regions  indi- 
cate there  were  several  stimuli  in  a given  time  period  when 
few  were  actually  presented.  These  results  give  rise  to  the 
conclusion  that  some  components  of  the  on-going  activity  match 
well  with  the  filter.  It  is  suggested  that  the  effectiveness 


v 


of  the  matched  filter  approach  may  be  improved  by  prefilter- 
ing the  sample  data  to  remove  the  contaminating  components. 
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IMPLEMENTATION  OF  A METHOD  TO 
DETECT  THE  SINGLE  VISUAL  EVOKED  RESPONSE 

I.  Introduction 

The  visually  evoked  response  (VER)  is  one  neural  response 
to  an  external  event.  In  practice,  discrete  visual  stimuli 
such  as  light  flashes  can  he  discerned  in  the  electroenceph- 
alogram (EEG),  and  these  representations  provide  valuable  in- 
formation concerning  the  processing  of  the  visual  stimulus. 
However,  the  neural  event  generated  in  response  to  this  stim- 
ulus is  extremely  small  with  respect  to  the  on-going,  spon- 
taneous electrical  activity  of  the  brain  as  measured  by  the 
EEG,  creating  a signal-to-noise  ratio  which  is  unfavorable. 

The  typical  way  to  improve  this  ratio  has  been  to  average  the 
response  to  a great  number  of  stimulus  presentations,  allowing 
the  average  of  the  random  activity  to  build  slowly,  while  the 
actual  evoked  response  builds  quickly.  While  satisfactory  in 
many  applications,  this  technique  precludes  using  the  evoked 
response  in  a number  of  other  applications.  In  particular, 
there  are  instances  in  which  the  VER  is  useful  only  if  it  can 
be  detected  after  one  stimulus  presentation. 

Purpose 

The  purpose  of  this  research  was  to  implement  a method  of 
detecting  the  visually  evoked  response  after  a single  stimu- 
lus presentation.  The  scope  of  this  work  is  limited  to  the 
signals  and  signal  detection  scheme  used  in  the  implementa- 
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tion.  However,  a brief  discussion  of  the  means  of  obtaining 
the  experimental  data  is  in  order  and  will  be  presented  in 
Chapter  II. 

Background 

Approach.  The  approach  taken  was  to  consider  signal 
detection  methods  known  from  information  processing  theory. 
The  matched  filter  technique  was  selected  since  response 
detection,  and  not  reproduction  (fidelity)  of  the  response 
was  the  primary  objective. 

A matched  filter  maximizes  the  ratio  of  the  peak  signal 
to  rms  noise  and  is  used  to  detect  the  presence  of  the  in- 
dividual evoked  response  following  a stimulus  presentation. 
With  the  matched  filter  the  output  signal-to-noise  ratio  is 
a function  solely  of  the  energy  in  the  signal  (VER)  and  the 
spectral  density  of  the  on-going  background  activity  which 
will  be  referred  to  as  noise.  It  is  the  energy  in  the  sig- 
nal that  provides  its  ultimate  detectability  in  noise  (Ref. 
iiM6). 

The  filter  was  constructed  (e.g.,  simulated  on  a digital 
computer)  from  a template  of  the  VER.  When  an  individual 
evoked  response  to  a single  stimulus  presentation  closely 
matches  the  template  the  filter  output  is  maximized.  That 
is,  the  filter  output  is  a pulse  which  indicates  that  the 
stimulus  was  perceived  by  the  brain.  In  essence,  the 
matched  filter  procedure  is  a cross-correlation  between  the 
expected  waveform  (or  template)  and  the  incoming  data.  When 
the  two  match  well,  the  resulting  cross-correlation  is  high. 
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If  a predetermined  threshold  is  exceeded  the  decision  is 
made  that  an  evoked  response  has  occurred. 

Assumptions.  The  implementation  of  the  approach  out- 
lined above  required  some  basic  assumptions  to  be  made. 

First,  it  was  assumed  that  the  signal  activity  will  be  simi- 
lar with  each  stimulus  presentation.  It  was  also  assumed 
that  the  averaged  VER  is  a valid  representation  of  each  in- 
dividual evoked  response.  A further  assumption  was  to  as- 
sume that  the  background  activity  and  evoked  response  are 
additive.  And,  finally,  the  noise  activity  and  the  VER  were 
assumed  to  be  uncorrelated.  The  validity  of  these  assump- 
tions will  be  discussed  along  with  the  results  in  Chapter  IV. 

Relationship  of  the  EEG  and  VER.  The  EEG  and  VER  are 
complex  phenomena  whose  specific  electrogeneses  are  unknown 
(Ref.  2i40).  The  background  activity  is  approximately  ten 
times  as  great  as  the  evoked  response  to  a single  visual 
stimulus  with  both  having  amplitudes  measured  in  microvolts. 

The  EEG  is  often  considered  to  be  a mixture  of  activity 
compartmentalized  into  various  frequency  ranges.  Starting 
with  the  low  frequency  range  there  is  the  delta  rhythm  or 
activity  which  is  usually  described  as  less  than  4 Hz.  Delta 
rhythm  does  not  usually  appear  in  normal  awake  adults.  Next, 
the  theta  activity  appears  in  the  frequency  band  of  4-7  Hz. 
This  activity  is  usually  of  low  voltage.  The  most  often  en- 
countered activity  is  alpha,  which  occurs  in  the  range  of  8- 
13  Hz.  This  activity  is  usually  large  and  quite  prominent. 
Alpha  or  alpha-like  activity  appears  to  be  the  major  contami- 
nant of  the  VER.  The  two  remaining  rhythms  are  beta, 


usually  described  as  low  voltage,  fast  BEG  in  the  range  of 
18-30  Hz,  and  gamma  which  occupies  the  frequency  range  of 
35-50  Hz.  This  latter  activity  is  seldom  encountered  or  de- 
scribed and  has  little  general  acceptance  (Ref.  3t25-26). 

The  data  acquisition  procedure  is  discussed  in  Chapter 
II.  A complete  discussion  of  the  signal  detection  metho- 
dology follows  in  Chapter  III.  Results  are  presented  and 
discussed  in  Chapter  IV.  An  outline  of  the  conclusions 
drawn  and  recommendations  for  future  work  are  offered  in 
Chapter  V. 


II.  Data  Acquisition 


Data  were  obtained  via  a VER  experiment  designed  by- 
Major  Robert  O'Donnell  of  the  Aerospace  Medical  Research 
Laboratory  at  Wright-Patterson  Air  Force  Base,  Ohio.  The 
experiment  was  designed  to  optimize  the  averaged  evoked  re- 
sponse and  thus  increase  the  chances  for  successful  detec- 
tion. The  experimental  procedures  and  instrumentation  are 
described  in  great  detail  by  a previous  Air  Force  Institute 
of  Technology  thesis  (Ref.  4- 1 8—33 ) • In  running  a VER  experi- 
ment there  are  many  variables  to  be  considered  which  influ- 
ence the  VER.  Perry  and  Childers  present  a detailed  discus- 
sion of  these  variables  and  their  effect  on  experimental 
results  (Ref.  3«^9-84). 

In  its  basic  form  the  VER  experiment  consists  of  the 
following  steps.  Electrodes  in  contact  with  the  scalp  con- 
duct the  continuous  electrical  activity  to  an  amplifier. 

The  amplified  signal  is  then  fed  to  a signal  averaging  com- 
puter. If  the  signal  were  recorded  prior  to  entering  the 
computer,  it  would  be  described  as  the  EEG.  In  order  to 
obtain  a VER,  the  computer  triggers  the  light  stimulus  some 
predetermined  number  of  times  and  simultaneously  averages 
the  scalp  electrical  activity  evoked  between  stimuli.  The 
resulting  average  of  the  nonrandom  activity  is  the  VER.  A 
control  or  noise  trial  is  obtained  in  exactly  the  same  man- 
ner, except  the  stimulus  is  occluded. 
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Data  Description 

For  purposes  of  this  research  four  channels  of  informa- 
tion were  recorded.  One  channel  consists  of  vertex  data  with 
a second  channel  containing  the  occipital  information.  The 
third  and  fourth  channels  contained  timing  and  stimulus  trig- 
gering data.  All  follow-on  data  processing  was  duplicated  to 
determine  which  region  produced  the  higher  probability  of 
evoked  response  detection. 

The  experiment  consisted  of  two  parts.  During  the  first 
series  the  template  was  constructed.  In  the  second  part  test 
data  were  recorded  which  would  allow  study  of  the  background 
activity  while  at  the  same  time  providing  sample  data.  This 
sample  data  would  later  be  used  to  test  the  matched  filter. 

Template  Construction.  Eight  averaged  VER's  were  ob- 
tained for  each  channel.  The  stimulus  rate  used  was  one  pre- 
sentation every  five  seconds  for  eight  periods  of  five 
minutes  each.  A black  and  white  checkerboard  pattern  was 
used  as  the  stimulus.  The  individual  responses  for  each 
period  were  computer  averaged  to  arrive  at  the  eight  aver- 
aged VER's  mentioned  previously.  One  template  for  each  chan- 
nel was  then  computed  by  taking  the  average  of  these  eight. 
The  resulting  templates  are  shown  in  Figs.  1 and  2 respec- 
tively. Each  template  is  approximately  20  microvolts  peak- 
to-peak  (maximum)  and  500  milliseconds  in  duration. 

Sample  Data.  The  sample  data  were  obtained  by  running 
three  10  minute  sessions  on  the  day  after  the  templates  were 
constructed.  A random  stimulus  presentation  rate  was  used 
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VERTEX  (CM.  1) 


TIME  CMS) 


TIME  CMS) 
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for  each  of  these  sessions.  For  the  first  run  the  rate  was 
an  average  of  one  stimulus  every  60  seconds.  Next,  an  aver- 
age rate  of  one  every  30  seconds  was  used,  and  for  the  final 
session  one  every  10  seconds.  Since  the  evoked  response  has 
little  effect  on  the  background  activity  the  combined  30 
minutes  of  sample  data  were  assumed  to  approximate  raw  EEG 
for  purposes  of  noise  computations.  Examples  of  vertex  and 
occipital  data  are  shown  in  Figs.  3 and  4 respectively.  The 
time  period  for  each  sample  is  six  seconds  and  the  approxi- 
mate peak-to-peak  maximum  amplitude  is  100  microvolts. 

Analog-to-Digltal  Conversion 

All  data  were  digitized  on  the  hybrid  computer  facility 
of  the  Computer  Science  Center  at  Wright-Patterson  Air  Force 
Base,  Ohio.  Personnel  of  the  Analog  and  Hybrid  Branch  devel- 
oped programs  to  digitize  the  data,  convert  it  to  a 60-bit, 
CDC  6600  compatible  representation,  and  store  it  on  magnetic 
tapes.  A sampling  rate  of  1000  per  second  was  used  which 
far  exceeds  the  minimum  requirement  to  sample  at  twice  the 
highest  signal  frequency. 
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III.  Matched  Filter  Design  and  Implementation 


The  fundamental  problem  to  solved  is  the  detection  of 
a known  signal  that  is  buried  in  noise.  Its  solution  calls 
for  a filter  which  will  maximize  the  signal-to-noise  ratio 
at  its  output.  Under  the  assumptions  listed  in  Chapter  I 
matched  filtering  is  the  optimum  method  of  implementing  this 
solution. 

This  chapter  begins  with  an  overview  of  the  matched  fil- 
ter derivation.  Following  this  is  a discussion  of  the 
matched  filter  simulation  using  the  digital  computer.  In- 
cluded in  this  discussion  are  a description  of  the  discrete 
Fourier  transform  (DFT)  and  the  method  used  to  determine  the 
noise  power  density  spectrum.  And,  finally,  the  filtering 
algorithm  used  to  filter  the  sample  data  is  presented. 

Matched  Filter  Derivation 

Considering  the  VER  template  as  the  signal  s(t),  its 
frequency  spectrum  S(f)  is  given  by  the  Fourier  transform 

S(f)  = 7 s(t)e_j2lTftdt  . (1) 

— OO 

The  signal  is  applied  with  additive  noise,  n(t)  (the  on- 
going background  activity),  to  the  input  of  a linear  filter 
with  transfer  function  H(f).  It  is  desired  to  determine  the 
characteristics  of  H(f)  that  will  give  the  greatest  ratio  of 
peak  signal  to  rms  noise  at  the  output  when  the  input,  x(t), 
is  s(t)  + n(t).  The  linear  aspects  of  the  filter  allow  the 
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sigral  and  noise  to  be  treated  separately  and  the  output 
results  added. 


V 


If  signal,  s(t),  alone  is  applied  to  the  filter,  the 
output  , g(t),  is  given  by 


g(t)  = 7 S(f )H(f )ej2nft  df  . (2) 

— 00 

The  magnitude  of  g(t)  will  be  maximum  at  some  time  t = tQ,  or 

I «(*„>  I = I 7 S(f)H(f)ej2Trfto  df  I > g(t)  . (3) 

If  the  input  noise  power  density  is  considered  to  be  B(f) 
watts/Hz,  the  mean  square  noise  power  at  the  output  of  the 
filter,  Nq,  is  seen  to  be 

Nq  = 7 B(f ) | H(f)  |2  df  (4) 


where , 


B (f ) = | A(f ) |2  . (5) 

The  ratio  to  be  maximized  by  the  proper  selection  of  H(f)  is 
the  peak  signal  to  rms  output  noise,  | g(tQ)  | /VNq  . This 
is  accomplished  by  squaring,  dividing  by  the  energy  in  the 
input  signal,  E,  and  applying  Schwartz’s  inequality.  The 
steps  required  are  shown  below. 
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A(f ) 


I g(t  ) I2  I s S(f)H(f)ej2Trfto  df  |; 

' 01  ' -oo  A(f ) ' 


_ oo  _ 

s ' (f)  | H (f ) |2  df  S | S(f ) |2  df 

mm  Y)  mm  OO 

oo  S (f ) ipirft  1 2 

J*  H(f )A(f )ej2"fto  df  2 

1 A(f ) L 

S B(f ) | H(f ) |2  df  J*  | S (f ) | 2 df 


J I — ^ I df  J*  I H(f  )A(f  )ej2Trfto  I2  df 
1 A(f)  1 --  1 1 


< — 
— oo 


S B (f ) | H(f ) |2  df  S | S(f ) y df 


The  ratio  can  therefore  be  maximized  by  letting 


S (f) 

H(f  )A(f ) = -g , 

A (f) 


H (f ) = 


S*(f)e-J2yft< 

A(f)A*(f) 


S*(f)e“^ftc 
B(f ) 


The  filter  may  have  an  attenuation  or  gain,  K,  that  is  not  a 
function  of  frequency.  Thus,  for  a non-white  noise  power 


density  spectrum,  the  transfer  function  of  the  matched  filter 
is, 

KS*(f  )e"^2Tffto 

H(f)  * . (12) 

B ( f ) 

Equation  (12)  demonstrates  that  the  optimum  filter  weights 
most  heavily  those  spectral  components  of  the  input  signal 
which  greatly  exceed  the  noise  spectral  density  and  de- 
emphasizes  the  spectral  regions  in  which  the  noise  spectral 
density  is  much  larger  than  the  signal  spectrum  (Ref.  5«146). 

Matched  Filter  Simulation 

The  output  of  the  analog-to-digital  process  described  in 
Chapter  I is  a discrete  time  signal  represented  by  a digital 
sequence  x(n).  For  this  sequence  it  is  possible  to  utilize 
an  alternative  Fourier  representation  referred  to  as  the 
discrete  Fourier  transform  (DFT ) . The  DFT  is  a Fourier 
representation  of  a finite-length  sequence  which  is  itself  a 
sequence  rather  than  a continuous  function.  The  DFT  is  es- 
sential for  the  computer  implementation  of  a variety  of  digi- 
tal signal  processing  problems.  Computation  time  for  the 
DFT  has  been  reduced  significantly  by  using  fast  Fourier 
transform  (FFT)  algorithms  such  as  the  one  used  in  this 
study  (Ref.  6 « 331 ) • A detailed  discussion  of  the  DFT  and 
derivation  of  computer  algorithms  for  the  FFT  is  presented 
in  Reference  (6). 

The  discrete  Fourier  transform  of  a sampled  continuous 
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time  function  f(t)  is 


where, 


N-i 

F(kv)  = E f (nT)e“JvTnk 
n=0 


(13) 


f(nT)  are  samples  of  f(t),  0 < n < N-l 

T is  the  sampling  interval 

F(kv)  are  samples  of  F(f),  0 < k < N-l 

v = 2tt/NT,  and  N = number  of  samples  of  f(t). 
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Its  inverse  discrete  Fourier  transform  (IDFT)  is  given  by 

f(nT)  = ^ N£1F(kv)ejvTnk  . (14) 

w k=0 

Since  the  data  were  sampled  at  a rate  greater  than  twice  the 
highest  frequency  of  F(kv),  the  discrete-time  Fourier  trans- 
form of  f(nT)  is  equivalent  to  samples,  equally  spaced  in 
frequency,  of  the  Fourier  transform  of  the  sequence  f(t). 

The  discrete  matched  filter  transform  can  therefore  be 
written  as 

•*<»>  ■ *§7#-  <!5> 

where  B(w)  is  the  power  density  spectrum  of  the  background 
activity.  This  is  the  discrete  form  of  Eq.  (12)  with  the 
exponential  term  dropped.  The  exponential  serves  to  specify 
the  time  tQ  where  the  peak  output  will  occur.  However,  the 
delay  inherent  in  the  filter  can  be  observed  in  the 
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convolution  of  each  template  with  its  corresponding  matched 
filter.  For  both  cases  t = 512  milliseconds  as  can  be  seen 
in  Figs.  5 and  6. 

The  problem  at  this  point,  then,  is  to  compute  the  DFT  of 
the  filter.  In  order  to  do  this  the  background  power  density 
spectrum  must  first  be  derermined. 

Noise  Power  Spectrum  Estimation.  The  background  activity 
is  of  indefinite  duration  and  not  periodic.  Consequently, 
the  discrete  time  signal  which  represents  this  activity  can 
be  modeled  best  in  terms  of  an  inf inite-duration  infinite- 
energy  signal.  Many  of  the  properties  of  this  type  of  signal 
can  be  summarized  in  terms  of  a finite-energy  sequence  called 
the  autocovariance  sequence,  for  which  the  Fourier  transform 
exist.  As  will  be  developed  below,  the  determination  of  the 
autocovariance  sequence  and  its  Fourier  transform  leads  to  a 
useful  interpretation  in  terms  of  the  frequency  distribution 
of  the  power  in  the  signal.  For  an  exact  signal  the  power 
density  spectrum,  P(w),  is  defined  as  the  Fourier  transform  of 
its  autocovariance  sequence. 

The  fundamental  concept  in  the  mathematical  representa- 
tion of  infinite-energy  signals  is  that  of  a random  process. 

In  the  discussion  that  follows,  it  is  assumed  that  the  reader 
is  familiar  with  fundamental  concepts  of  probability  and 
random  processes. 

T*e  procedure  used  in  determining  a useful  estimate  of 
the  noise  power  spectrum  is  basically  the  Welch  method  (Ref. 

6i 553-55*0 • Implementation  of  this  method  calls  for  averaging 
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modified  periodograms . Periodograms  are  themselves  estimates 
of  the  spectrum  of  a finite  sequence.  The  usefulness  of  the 
noise  power  spectrum  so  determined  requires  that  it  converge 
to  some  meaningful  result  as  the  measurement  interval  in- 
creases. A simple  and  useful  way  of  defining  convergence  is 
by  requiring  that  the  variance  of  the  estimate  approach  zero 
as  the  measurement  interval  increases  (Ref.  7«399)» 

As  an  estimate  of  the  power  density  spectrum  consider  the 
Fourier  transform  of  the  biased  auto^ovariance  estimate  c(m), 
e.g. , 


IN  (w)  = V c(m)e"jwn  . (16) 

N m = — ( N— 1 ) 

Since  the  Fourier  transform  of  the  real  finite-length 
sequence  x(n),  0 < n < N-l , is 

X(ejw)  = Z1x(n)e"'jwn  (17) 

n=0 

and  it  can  be  shown  that 

IN  (w)  = - | X(ejw)  |2  . (18) 

The  spectrum  estimate  ( w ) is  called  the  periodogram  (Ref. 

6 « 542).  However,  the  variance  of  IN(w)  does  not  approach 
zero  as  N approaches  infinity. 
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var  pN(w)] 


c°v  pN(w),IN(w)J 
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- p2(*>  [*  +CaifSj}]2]  <«> 

Thus,  the  periodogram  is  not  a consistent  estimate. 

A standard  approach  to  reducing  the  variance  of  estimates 
is  to  average  over  a number  of  independent  estimates.  The 
application  of  this  approach  to  spectrum  estimation  is  often 
attributed  to  Bartlett  (Ref.  6«5^8).  Therefore,  the  data 
sequence  x(n),  representing  the  background  activity,  is  di- 
vided into  K segments  of  M samples  each  so  that  N=KMj  e.g., 
forming  the  segments 


.(i) 


(n)  = x(n+im-m),  0 < n < m-1 , 1 < i < K . (20) 


Next,  the  K periodograms  are  computed  from  Eq.  (21). 


(w)  = _ 
(vf)  M 


ML1x(i)(n)e-^wn 

n=0 


, 1 < i < K (21) 


The  spectrum  estimate  is  then  defined  as, 


1 K 


B(w)  - if  • 


(22) 


To  determine  if  this  is  a useful  estimate  consider  its 
variance  as  described  in  Eq.  (23)* 
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var  [b(w)^| 


= K Var  [XM(W)] 

+Crtt]2]  <23> 

It  is  clear  that  the  variance  of  B(w)  is  inversely  propor- 
tional to  the  number  of  periodograms  averaged,  and  as  K gets 
large,  the  variance  approaches  zero.  Therefore,  the  Bartlett 
estimate  is  a consistent  estimate. 

The  Welch  method  is  a modification  of  the  Bartlett  pro- 
cedure. It  is  particularly  well  suited  for  direct  computa- 
tion of  the  power  spectrum  estimate  using  the  FPT.  In  this 
procedure  the  data  record  is  again  sectioned  into  K = N/M  seg- 
ments of  M samples  each  as  defined  in  Eq.  (20).  However,  the 
periodograms  are  modified  by  applying  the  triangular  window 
sequence,  w(n),  directly  to  the  data  segments  before  computa- 
tion of  the  periodogram.  Thus,  the  K modified  periodograms 
are  defined  by 


1 

MU 


ME1x(i)(n)w(n)e"jwn 

n=0 


i-1.2. 


>K 


(24) 


where , 


w(n) 


2n 

N-l  ’ 


0 < n < 


N-l  . . 

— < n < 


N- 


1 
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and, 


1 M’1  2 

U = i E *r(n)  . (25) 

n=0 

The  spectrum  estimate  is  then  defined  as, 

B(w)  = £ I ^(w)  . (26) 

Welch  (Ref.  8)  shows  that  if  the  segments  of  x(n)  are  non- 
overlapping, such  as  in  the  present  case,  then 

var  Qb  (w )^J  « ^ P2(w)  . (27) 

This  is  the  same  as  for  the  Bartlett  method  and  is  therefore 
a consistent  estimate. 

As  described  above,  then,  a consistent  estimate  of  the 
noise  power  spectrum  was  computed  by  averaging  modified 
periodograms.  The  sample  data  were  divided  into  3510  seg- 
ments of  512  samples  each.  Each  segment  was  multiplied  by  the 
512  point  window  sequence  and  then  a 512  point  PFT  algorithm 
was  computed.  The  magnitude  of  the  FFT  for  each  segment  was 
squared  and  the  power  spectrum  P(w)  was  determined  by  taking 
the  average  of  these  squares.  Figures  7 and  8 show  the  spec- 
trum for  the  vertex  and  occipital  channels  respectively.  The 
magnitude  shown  is  in  logj0  format. 

Filter  Response.  The  matched  filter  response  was  com- 
puted from  Eq.  (15)  by  first  computing  H(w).  Dividing  the 
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complex  conjugate  of  S(w)  by  B(w)  yields  H(w).  The  constant 
K was  arbitrarily  set  equal  to  one.  Thus,  the  filter  re- 
sponse, h(n),  is  computed  by  finding  the  IDFT  of  H(w). 

The  first  computation  of  h(n)  did  not  produce  the  desired 
results.  Upon  inspection  of  Figs.  7 and  8,  the  higher  fre- 
quency components  are  seen  to  be  over  emphasized.  To  correct 
this  condition  the  power  spectrum  for  both  channels  were 
modified.  The  curves  were  first  mathematically  smoothed. 
Then,  to  reduce  the  peak  magnitudes  the  eighth  root  of  the 
power  spectrum  was  computed.  The  results  of  this  procedure 
are  shown  in  Figs.  9 and  10. 

The  matched  filter  response  for  each  channel  was  then 
recomputed.  As  shown  in  Figs.  11  and  12,  the  resulting  re- 
sponse curves  very  closely  approximate  the  template  signals 
running  backwards  in  time. 

Filter  Implementation 

In  the  general  case  the  filtered  output,  y(n),  is  com- 
puted by  performing  a linear  convolution  of  the  input  signal 
and  filter  response,  e.g., 

N-l 

y(n)  = £ x(m)h(n-m)  (28) 

n=0 


where , 


x(n)  is  the  input  signal 
h(n)  is  the  filter  response. 
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If  both  x(n)  and  h(n)  are  N-point  sequences  the  output 
sequence  will  have  at  most  2N-1  nonzero  points. 


For  the  problem  at  hand  the  input  sequence  is  the  on-going 
EEG  containing  randomly  occurring  VER  signals.  This  sequence 
has  an  indefinite  length  and  therefore  must  be  segmented.  The 
segment  size  selected  was  512  points  which  is  equal  in  length 
to  the  filter  response  sequence. 

The  discrete  Fourier  transform  was  utilized  to  effect  the 
required  linear  convolution  by  performing  a circular  convolu- 
tion between  each  input  data  segment  and  the  filter  response. 

A circular  convolution  is  implemented  by  computing  the  DFT's 
of  the  two  sequences,  multiplying,  and  computing  the  IDFT  of 
the  product.  To  effect  the  linear  convolution,  however,  there 
are  some  additional  steps  involved. 

Prior  to  computing  the  DFT's  the  sequences  must  be  zero 
filled  to  at  least  2N-1  noints.  In  this  case  each  sequence 
was  made  1024  points  in  length  and  the  DFT's  were  then  com- 
puted on  this  basis.  To  construct  the  filtered  output  for  the 
entire  input  sequence  the  results  for  the  individual  filtered 
segments  must  be  added  in  a particular  manner. 

The  procedure  used  to  construct  the  filtered  output  from 
filtered  segments  is  referred  to  as  the  overlap-add  method 
(Ref.  6ill2-113).  Each  filtered  segment  of  length  1024 
noints  will  ovr  ap  the  previous  segment  by  512  points.  Thus, 
each  filtered  segment  is  overlapped  with  and  added  to  the  pre- 
vious segment  to  construct  the  filtered  output.  The  program 
used  to  implement  this  procedure  is  included  in  Appendix  A. 
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Computer  memory  space  requirements  were  high.  To  filter 
30  seconds  of  data  required  approximately  200,000  octal  words 
and  consequently  this  was  the  maximum  amount  of  data  filtered 
per  run.  Results  of  the  filter  implementation  are  presented 
in  the  following  chapter. 
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IV.  Results  and  Discussion 


Using  the  computed  filter  response  and  the  overlap-add 
filtering  scheme  described  in  Chapter  III,  30  second  segments 
of  sample  data  were  filtered  for  both  the  vertex  and  occipital 
channels.  The  filtered  output  should  show  a positive  peak 
whenever  a stimulus  was  perceived  by  the  subject.  This  peak 
should  be  well  defined  and  have  greater  amplitude  than  the 
output  where  no  stimulus  occurred. 

The  sample  data  were  arranged  in  three  10  minute  files 
according  to  the  stimulus  rate  used.  File  one  is  the  one  per 
60  seconds  ratej  file  two  contains  one  stimulus  per  30  seconds} 
and  file  three  contains  the  one  per  10  seconds  rate. 

Data  from  each  file  were  filtered  producing  similar  results 
in  all  cases.  Therefore,  the  discussion  that  follows  will  be 
limited  to  the  results  of  filtering  the  first  90  seconds  from 
file  one  for  both  channels.  Because  the  stimulus  presenta- 
tion rate  is  an  average  rate  there  are  six  stimuli  occurring 
in  these  90  seconds.  These  occurred  at  0,15,30,^5,60,  and 
76.5  seconds  into  the  experiment.  To  simplify  the  discussion 
these  will  henceforth  be  referred  to  as  stimulus  A,  B,  C,  D, 

E,  and  F respectively. 

As  will  become  obvious  upon  examination  of  the  filtered 
output  data,  a decision  as  to  whether  or  not  a given  evoked 
response  was  detected  is  somewhat  arbitrary.  The  criteria 
for  judging  a positive  detection  was  based  primarily  upon 
the  appearance  of  a sharp  peak  .5  to  1.0  seconds  after  the 
stimulus  was  presented.  First  to  be  discussed  is  the  vertex 
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data  results  followed  by  the  results  for  the  occipital  data. 


Results  for  Vertex  Data 

The  vertex  filter  output  is  shown  in  Figs.  13A,  B,  and  C. 
Clearly,  the  desired  pattern  is  not  in  evidence.  There  are 
some  peaks  which  are  distinct  and  of  greater  amplitude  than 
a majority  of  the  remaining  output.  However,  these  do  not 
coincide  with  the  presentation  of  stimuli.  Closer  examination 
of  the  smaller  peaks  reveal  concurrence  (within  .5  to  1.0 
seconds  afterward)  with  B,  C,  D,  E,  and  F.  The  amplitude  of 
these  peaks  ranges  from  2.k  to  5*0  microvolts.  Although  the 
match  between  the  filter  and  some  of  these  single  VER's  is 
relatively  poor  they  do  indicate  a positive  detection  rate  of 
83.3  %•  This  is  speculative  to  be  sure.  It  is  apparent  that 
the  on-going  background  activity  contains  other  components 
that  match  well  with  the  matched  filter.  For  instance,  the 
high  amplitude  peaks  appearing  at  approximately  39  and  80 
seconds  coincide  with  alpha  or  alpha-like  rhythm  in  the  back- 
ground activity. 

Results  for  Occipital  Data 

The  occipital  filter  output  is  shown  in  Figs.  \kk,  B,  and 
C.  As  with  the  vertex  filter  the  desired  output  was  not 
achieved.  Again  there  are  some  pe.- ..s  which  coincide  with  the 
presentation  of  stimuli.  Specifically,  these  appear  for  A,  B, 
C,  D,  and  F with  amplitudes  ranging  from  1.75  to  5*0  micro- 
volts. This  would  also  indicate  a positive  detection  rate  of 
83.3  %•  It  should  be  noted  that  high  amplitude  peaks  occur  at 
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39  and  80  seconds  as  in  the  vertex  output.  However,  there  are 
other  high  amplitude  peaks,  indicating  a good  match,  such  as 
the  one  at  approximately  11  seconds  for  which  no  apparent  ex- 
planation is  available. 
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V.  Conclusions  and  Recommendations 


Conclusions 


The  objective  of  this  project  was  to  implement  a detection 
scheme  which  would  detect  the  visually  evoked  response  to  a 


single  stimulus.  This  VER  is  known  to  be  buried  in  the  brain's 


on-going  electrical  activity.  Based  on  the  results  of  this 


study  the  following  conclusions  were  drawn i 


1.  The  matched  filter  will  yield  a positive  detec- 
tion rate  of  over  80  %.  However,  its  usefulness  is  nullified 


by  the  high  number  of  false  alarm  indications  in  the  filter 


output. 


2.  At  this  point  no  distinction  can  be  made  as  to 


whether  the  vertex  data  or  the  occipital  data  is  more  prefer- 


rable  for  single  VER  detection. 


3.  Because  of  the  additional  apparent  matchings  the 


matched  filter  by  itself  is  not  an  effective  method  of  single 


evoked  response  detection. 


Recommendations 


It  is  recommended  that  further  study  in  this  area  consider 


some  modifications  to  the  basic  matched  filter  approach  de- 


scribed in  this  report.  The  validity  of  the  averaged  VER, 


which  is  used  as  the  template,  requires  further  study.  There 


can  be  marked  differences  between  any  two  individual  responses 


although  stability  of  the  final  averaged  VER  does  not  suggest 


the  fluctuation  of  these  contributing  responses  (Ref.  3.11-12, 


2.73). 


kZ 


Further  study  is  required  into  the  relationship  between 
the  VER  and  the  on-going  background  activity.  The  inter- 
action between  them  is  additive  only  to  a first  approximation. 
There  is  evidence  of  a non-linear  or  non-stationary  inter- 
action between  the  background  activity  and  an  evoked  response. 
This  non-linear  interaction  can  also  be  inferred  from  an  ex- 
amination of  evoked  responses  to  single  flashes  of  light  such 
as  in  the  present  case.  A large  percentage  of  the  time  these 
individual  evoked  responses  can  be  observed  by  the  eye,  de- 
spite the  level  of  background  activity  - that  is,  the  noise 
seems  to  be  suppressed  during  the  major  portion  of  the  evoked 
response  (Ref.  3*27-28). 

A final  recommendation  is  to  study  the  effectiveness  of 
prefiltering  the  sample  data  with  a bandpass  filter  to  elimi- 
nate the  low  and  high  frequency  noise  components.  Also,  since 
alpha  or  alpha-like  activity  appears  to  be  the  major  contami- 
nant of  the  VER  more  than  one  bandpass  filter  may  be  in  order. 
There  is  some  speculation  that  comb  filtering  may  provide  some 
positive  results.  This  technique  utilizes  a series  of  bandpass 
filters  which  pass  only  those  frequency  components  that  contri- 
bute to  the  evoked  response  and  reject  those  frequency  re- 
gions which  are  conjectured  to  be  contaminated  by  noise  alone 
(Ref.  3*113). 
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visual  evoked  response.  The  approach  taken  is  based  on  the 
similarity  of  the  problem  to  that  of  radar  detection  in 
which  a known  signal  is  buried  in  noise.  Using  the  discrete 
Fourier  transform  (DFT)  and  other  digital  signal  processing 
techniques  a filter  matched  to  the  VER  template  was  computed 
and  implemented  on  a digital  computer. 

Computation  of  the  matched  filter  required  an  estimate 
of  the  noise  power  spectrum,  where  the  noise  in  this  case  is 
the  on-going  electrical  activity  of  the  brain.  This  was  ac- 
complished by  averaging  modified  periodograms . Data  from 
both  the  vertex  and  occipital  regions  were  filtered  to  deter- 
mine if  one  region  produced  better  results  than  the  other. 

The  results  indicate  that  the  matched  filter  will  detect 
the  single  evoked  response.  However,  the  false  alarm  rate 
is  high.  The  filtered  output  for  data  from  both  regions 
indicate  there  were  several  stimuli  in  a given  time  period 
when  few  were  actually  presented.  These  results  give  rise 
to  the  conclusion  that  some  components  of  the  on-going  activ- 
ity match  well  with  the  filter.  It  is  suggested  that  the 
effectiveness  of  the  matched  filter  approach  may  be  improved 
by  prefiltering  the  sample  data  to  remove  the  contaminating 
components. 
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